Comment on "First-order phase transitions: equivalence between bimodalities 

and the Yang-Lee theorem" 
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I discuss the validity of a result put forward recently by Chomaz and Gulminelli [Physica A 
330 (2003) 451] concerning the equivalence of two definitions of first-order phase transitions. I show 
that distributions of zeros of the partition function fulfilling the conditions of the Yang-Lee Theorem 
are not necessarily associated with nonconcave microcanonical entropy functions or, equivalently, 
with canonical distributions of the mean energy having a bimodal shape, as claimed by Chomaz and 
Gulminelli. In fact, such distributions of zeros can also be associated with concave entropy functions 
and unimodal canonical distributions having affine parts. A simple example is worked out in detail 
to illustrate this subtlety. 
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Chomaz and Gulminelli Q have studied recently the 
equivalence of two different definitions of first-order phase 
transitions — one based on the nonconcavity of the micro- 
canonical entropy function and another based on the dis- 
tribution of the zeros of the partition function. My goal 
here is to question and correct one of their results with 
the help of a counterexample which I will then explain 
using some basic results of convex analysis. To start, let 
me summarize the main results found in Ref . 0] . The no- 
tation I will be using throughout is less general than the 
one used in it is simpler, but captures nevertheless 
the essence of the problem. 

Consider an n-body system with energy U and mean 
energy u = U/n. The partition function of the system is 
defined as 



Z n (f3) = / n n (u)e 



-{3nu 



(i) 



where VL n {u) represents the density of microstates with 
mean energy u. It is well-known from the work of Yang 
and Lee [2| that one way to make sense of nonanalytic 
points of the canonical free energy function 
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m £„(/?), 



(2) 



which signal the onset of phase transitions, is to study 
the distribution of complex zeros of Z n (Jf) in the limit 
n — ► co. In the case of first-order phase transitions, in 
particular, it is known that, as n — > oo, there is an ac- 
cumulation of zeros of Z n (J$) around some positive real 
value P c of the inverse temperature corresponding to the 
value at which <p{(3) is nondifferentiable, and that the loci 
of zeros in the vicinity of j3 c is parallel to the imaginary 
axis. This phenomenon is what Chomaz and Gulminelli 
refer to as the Yang-Lee Theorem. I shall refer to it my- 
self as the Yang-Lee Condition (YLC). Thus we say that 
a first-order phase transition appears in the thermody- 
namic limit of the canonical ensemble when the zeros of 
Zjfi) satisfy YLC. 



Now, what Chomaz and Gulminelli purported to show 
in |]| is that YLC is equivalent to another definition of 
first-order phase transitions based on the bimodal shape 
of the density function fi n (it) (see [l| and references 
therein). They showed first that if f2 n (u) has a bimodal 
shape which persists as n — > oo, a condition which I shall 
refer to as the Bimodal Condition (BC), then YLC is sat- 
isfied. Actually, the quantity which they focused on was 
not fl n (u) but the canonical distribution p n p(u), given 
by 



Pn,f}(u) 



n n (u)e-P nu e -n[Pu-s(u)] 



Zn((3) 



Zn(P) 



(3) 



Yet since a bimodality of f2„ (u) translates into a bimodal- 
ity of p n ,^{u), and vice versa, it does not matter which 
quantity we refer to, and, for the purpose of the presen- 
tation, I shall stick to Q n (u). 

It should be noted in passing that the result "BC im- 
plies YLC" had already been proved by Lee [3( ■ The real 
novelty of |1J is to attempt to prove the converse result, 
namely that YLC implies BC. It is this second result that 
Chomaz and Gulminelli claim to have proven, but which 
cannot be true in fact, as exemplified by the following 
counterexample. 

Consider a density of microstates of the form 



fi„(u) 



lA u e [o, A] 

elsewhere, 



(4) 



with A > 0. The partition function for this form of O n (u) 
is trivially evaluated and has for expression: 



nA 

Z n ((3) = — 
n 
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(5) 



Setting Z n (f3) — and solving for j3 £ C, we find that 
the zeros of the partition function must solve the equation 



-/3nA 



to e 



-0nA 



1 with the exclusion of j3 



1,2, 



zeros of Z n (f3) to be given by (3k 



0. This is equivalent 
and so we find the 
= ±2nik/(nA), k = 
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FIG. 1: A generic nonconcave entropy function s(u) and its concave envelope s**(u), which is affine (with slope f3 c ) where s(u) 
is nonconcave. Both entropies have the same Legendre-Fenchel transform which corresponds to the free energy function (p{j3). 



1,2,.... In terms of the fugacity z 
re-expressed as 



/3 , these can be 
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Now comes the contradiction: £l n (u) does not satisfy 
BC, but the zeros of the partition do satisfy YLC. In 
fact, the zeros of Z n (f3) are all aligned on the imaginary 
axis and pinch the real axis at the critical value [3 C — 
as n — > oo. Following the Yang-Lee theory 0], we then 
know that Z n {(3) must develop a nonanalytic point at f3 c 
as n — > oo, which translates into a nondifferentiable point 
of tp(f3). This can be verified by a direct calculation of 
the free energy: 
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(7) 



The left- and right-derivatives of the free energy at (3 C = 
being equal to (p'(j3 = 0") = A, <p'([3 = 0+) = 0, 
respectively, we also find that the latent heat for this 
example is equal to A. 

The important conclusion that we reach from this sim- 
ple counterexample is clear: YLC does not imply BC in 
general, as claimed in . To determine what the correct 
implication should be, I shall recall at this point three im- 
portant results of equilibrium statistical mechanics and 
convex analysis (see |j,|5||): 

(i) If £l n (u) has a bimodal shape that persists when 
n — > oo, then s(u) must be nonconcave over some range 
of mean energy. The converse statement is also true. 

(ii) The free energy function tp(f3) is always the 
Legendre-Fenchel transform of the microcanonical en- 
tropy function s(u); in symbols, 



tp(P) — inf{/3u — s(u)}. 



(8) 



This holds no matter what the shape of s(u) is, be it 
concave or not. 

(iii) Regions of mean energies over which s(u) is non- 
concave or is affine (i.e., is a line) are indicated at the 
level of the canonical free energy <p(/3) by the existence 



of points of <p(/3) where this function is nondifferentiable. 
The fact that affine parts of s(u) also translate into non- 
differentiable points of can be understood by noting 
that s(u) and its concave envelope s**(u), defined by 



s**(u)=M{Pu-<p(0)}, 



(9) 



have the same Legendre-Fenchel transform, namely, 



<p(/3) = M{[3u - s(u)} = M{[3u - s**(u)}; (10) 

u u 

see Fig. ^ 

These results indicate altogether that a first-order 
phase transition in the canonical ensemble emerges from 
the point of view of the microcanonical ensemble in basi- 
cally two ways: either s(u) is nonconcave somewhere or 
else s(u) is affine somewhere |||. The first possibility is 
what was considered in 0, whereas the second possibil- 
ity is what I have considered in the counterexample, and 
what is precisely absent from . Hence at this point we 
have all the ingredients to conjecture what the correct 
relationship between YLC and BC is, namely: if the ze- 
ros of Z n (/3) satisfy LYC, then either fi n (it) is bimodal 
or else it has an affine part. Let me refer to the second 
possibility as the Affine Part Condition (APC). Then, 
the result is simply: YLC implies BC or APC. Combin- 
ing this result with what we already had, namely that 
BC implies YLC, we arrive finally at the following: YLC 
is satisfied if and only if BC or APC is satisfied. (This 
result should be understood as a conjecture rather than 
a rigorously proved result because, at this point, it re- 
mains to rigorously prove that YLC is a necessary and 
sufficient condition for <p{(3) to have a nondifferentiable 
point.) 

In the end, one may wonder whether entropies having 
affine parts are anything to worry about. There seems 
indeed to be a lack of structural stability inherent with 
such entropies, and so it is questionable whether they can 
show up in realistic models. This concern, however, is out 
of the scope of the present comment. Affine entropies 
present one with a theoretical possibility that one has 
to take into account when deriving general theoretical 
results. 
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